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The frustrated spin- 1/2 model with weakly anisotropic ferromagnetic nearest-neighbor and antifer- 
romagnetic next-nearest-neighbor exchanges is studied with use of variational mean-field approach, 
scaling estimates of the infrared divergencies in the perturbation theory and finite-size calculations. 
The ground state phase diagram of this model contains three phases: the ferromagnetic phase, 
- - - , the commensurate spin-liquid phase and the incommensurate phase. The non-trivial behavior of 

the boundaries between these phases and the character of the phase transitions in case of weak 
' anisotropy are determined. 

(N ; 

^ , I. INTRODUCTION 

The quantum spin chains with nearest-neighbor (NN) Ji and next-nearest-neighbor (NNN) interactions J2 have 
been a subject of numerous studies [T]. The model with both antiferromagnetic interactions Ji, J2 > (AF-AF model) 
I— — 1| is well studied ^, ^, Ji, ^ 6, 7]. Lately, there has been considerable interest in the study of F-AF model with the 
^ • ferromagnetic NN and the antiferromagnetic NNN interactions (Ji < 0, J2 > 0) [1, H, [13, 111 [S [H, H • One of the 
reasons is understanding of intri guing magnetic prop erties of a novel class of edge-sharing copper oxides which are 
described by the F-AF model [3, Il6l. ITtI. [isl . [l9l . [20| . In particular, these copper oxides show at low temperature a 
tendency to the formation of the incommensurate state with helicoidal magnetic ordering. 
, The Hamiltonian of the F-AF model is 
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where we put Ji = — 1 and J2 = J > and the periodic boundary conditions are imp lied. 

The isotropic case of this model (Ai = A2 = 1) is intensively studied last years [ill, [l^. [21I [2^. 12^ . It is known that 
the ground state of the isotropic version of the model (H]) is ferromagnetic at < J < 1/4 and it becomes a singlet 
incommensurate state for J > 1/4 [HjUHl- The phase transition at J = 1/4 is the second order one. 

The model with the anisotropy of exchange interactions is less studied, especially for the case of the small anisotropy. 
For example, the phase diagram of the model ([T]) with Ai = A2 has been studied in Ref.jl^ using the method of 
level spectroscopy. Unfortunately, this method becomes unreliable for J « 1/4 and Ai = A2 ~ 1 because of strong 
finite-size effects. 

, In real chain compounds the exchange interactions are anisotropic. The microscopic origin of these interactions is 
' the spin-orbit coupling. The indication on the anisotropy is a dependence of the saturation field on the direction of the 
t — \ external magnetic field [l3|. Though, as a rule, the anisotropy is weak (for example, for edge-shared cuprate LiCuV04 
' ESR detected a 6% anisotropy [27]), it can change the transition point from commensurate to incommensurate states 
^ ^ as well as the behavior of the model ([T]) in the vicinity of the transition point. Besides, the frustration parameter 
\J2/Ji\ — J estimated for some edge-sharing copper oxides is close to the quantum critical point 1/4 (for example, 
J ^ 0.28 — 0.3 for compound Li2ZrCuOi "28]). Therefore, taking into account both the frustration effects and the 
small exchange anisotropy near the transition point can be important for the analysis of the experimental data related 
to these compounds. 

In the isotropic case of ([T]) the ferromagnetic state is (TV -I- l)-fold degenerated at < J < 1/4. Weak easy-plane 
anisotropy Ai, A2 < 1 lifts this degeneracy and the ground state is in the sector with total 5*^ = at small J . One 
can expect that the increase of J induces the phase transition at some Jc to the incommensurate phase with 5*^ = 0. 
Besides, the character of this transition can be different from that in the isotropic case. 

In our analysis we focus on the behavior of the F-AF model ^ near the transition point from the commensurate to 
the incommensurate ground state and the influence of the weak anisotropic interaction on the T = phase diagram. 
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For simplicity we concentrate our attention on the particular case of the Hamiltonian ([T]) with A2 = 1 

H=- Y.{S:s:^, + SlSl^, + hS^Sl^, -\) + jY,{S^ ■ S„+2 - J) (2) 

(We added here constants for convenience.) 

However, we will show that the results for the model ^ with both Ai 7^ 1 and A2 ^ 1 are qualitatively similar to 
those for the model 

The paper is organized as follows. In Sec. II we consider a qualitative physical picture of the ground state phase 
diagram of the model ([2]) based on the classical approximation. In Sec. Ill we study the phase diagram of the model 
^ using the variational mean-field approach. The scaling estimates of the perturbation theory for the easy-plane case 
of the model ([2]) at J < 1/4 are presented in Sec. IV. In Sec.V we estimate infrared divergencies in the perturbation 
theory near the transition point J = 1/4. Sec. VI is devoted to the phase transition in the easy-axis case of the model 
([2]). In Sec. VII we present the phase diagram of the model ^ in the case Ai = A2 and summarize our results. 



II. CLASSICAL APPROXIMATION 



Let us start from the classical picture of the ground state of the model In the classical approximation the spins 
are vectors which form the spiral structure with a pitch angle ip between neighboring spins and canted angle 9 

S'^ — - cos((pn) sin6' 
— sa\{ipn) sin 9 

= ^cos0 (3) 



The classical energy per site is 

Eci{(p,9) _ 1 - A sin^l 



N 



[A - cos ^ - J(l - cos(2(^))] (4) 



The minimization of the energy ([4]) over the angles tp and 9 shows that there are three regions in (J, A) having 
different classical energies. In the region I ( J < i, A < 1) the energy is minimized by the choice of the angles tp — 
and 9 — ^. These angles correspond to the spin configuration with all spins pointing along the a;-axis and the energy 
is 

Ecl,^y = (5) 

In the region II (J < i, A > 1) and ( J > ;|, A — 1 > j ( J — i)^) the minimum of the energy is given by the angle 
9 — (and arbitrary ip). This is the fully polarized state with all spins up (or down) and the energy 

A-1 

EcL, = -N^— (6) 

In the region III ( J > i, A — 1 < -j ( J — i)^) the classical approximation shows helical spin structure in the x-y 
plane. The corresponding angles are 

-1 1 

ip = cos 

^ 4J 

« = i (7) 



and the classical ground state energy is 



N / ^ 1^ ^ 



The phase boundaries in the classical approximation for the model ^ are shown in Fig.l by thin dashed lines. 
One can see from Fig.l that the transition between phases I and II takes place on the isotropic line A = 1. This 
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FIG. 1: The phase diagram of the model (|2]). 

transition is a simple spin-flop, which is certainly of the first-order type. In the easy-axis case A > 1 the increase 
of the NNN exchange J leads to the first-order transition to the helical phase III with finite value of (/? (O, which 
is If = (8(A — l))^^"* at A ^ 1. In contrast to the easy-axis case, in the easy-plane part of the phase diagram the 
transition to the helical phase occurs at J — 1/4, where the pitch angle (fi = 0, indicating the second-order type of 
this transition. 



III. MEAN FIELD APPROACH 

To study the model ^ we use the variational mean-field approach (MFA) developed in Ref. [U, [2^ . According to 
this approach, we follow the classical picture and transform the local axes on n-th site by a rotation about the Z axis 
by <pn and then by a rotation about the Y axis by 9. The transformation to new spin-i operators 7y„ has a form 

S„ = R,i^n)Ry{0)r,n (9) 

where Ry{0) and Rz{ipn) are the operators of the corresponding rotations. 

The second step is the Jordan- Wigner transformation to the obtained Hamiltonian in terms of the rj operators. 
This transformation maps the Ty-spin model onto the model of interacting spinless fermions, which is then treated by 
the mean-field approximation including superconductor like correlations. The pitch and canted angles ip and are 
variational parameters in this approach. We omit here the details of this approach, because it is simple modification 
of that done in Ref. [2ll. |29| . and we present here only the results of this approximation. 

Generally, the phase diagram of the model ([2]) in the MFA contains the same three phases as predicted by the 
classical approximation and the boundaries between the phases are shown in Fig.l by thick dashed lines. In the 
region I (see Fig.l) the MFA shows the non-zero magnetization in the easy x — y plane. In the region II the fully 
polarized state |tt ■ ■ ■ T) represents the ground state. In the region III the MFA shows helical spin structure in the 
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X — y plane. However, as can be seen in Fig.l, the boundary between the phases I and III is substantially shifted. In 
the MFA this boundary in the vicinity of the point (J = A = 1) is approximately given by 

a « 8.0571-2^ (10) 

where a = 1 — A and 7 = J — ^ . 

The boundary between the phases II and III in the MFA is described by the equation 

Awl + 6.37^-'^ (11) 

Certainly, there is no LRO in the x-y plane in the phases I and III and in this respect the MFA is incorrect. 
However, the MFA gives a good estimate for the ground state energy in those phases. For example, in the phase I at 
J = the MFA reproduces correctly non-trivial critical exponent for the ground state energy 

5Eq w -0.0637Va^/2 (12) 
This estimate differs on 16% in numerical factor from the exact result [soj 

The MFA shows that the critical exponent 3/2 for the ground state energy remains up to the point J = 1/4, where 
the behavior of the ground state energy is changed to 

SEo ~ -0.07Na^/'^ (14) 

As was shown in Ref. fjlj the MFA gives also a good estimate for a critical exponent of the ground state energy in 
the isotropic case A = 1 of the helical phase III: 

SEo ~ -I.5857V7I2/7 (--^5) 

As will be shown below, the estimates of the ground state energies and the phase boundaries in the MFA given by 
Eqs.pO |l . ([TT|) are in a good accordance with scaling estimates and finite-size calculations. 



IV. PERTURBATION THEORY FOR EASY-PLANE CASE AT J < 1/4 

We are interested in the behavior of the model ([2]) in the vicinity of the isotropic case A = 1. For this aim it is 
natural to develop the perturbation theory (FT) in small parameter a = 1 — A 

H = Ho + Vj + Vs 

Ho = - ' ^"+1 ^ 4) 

Vj = j5](S„.S„+2-i) 

= a ^5,15,1+1 (16) 

At first let us consider the most simple case J = 0, where the ground-state energy at a ^ 1 is given by Eq. p^ . 
The ground state of Hq is ferromagnetic and is degenerated with respect to total . The perturbation Vq splits this 
degeneracy and in the first order in a we have 

{i,{Sn\V^msn) = ^^-^^^ (17) 

Thus, the first order shows that one should develop the FT for the lowest state \ipo) having total spin 5 = y and 
S"^ = 0. The perturbation series for the ground state energy can be written in a form: 

Eoia) = (^ol V^ + Vo. \ y„ + . . . IV'o) (18) 
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Suppose that the main contributions to the energy are given by low-lying excitations, which for isotropic ferromagnet 
with the spectrum = behave as 

Ek -Eo^ N-^ (19) 

The higher orders of the perturbation series contain more dangerous denominators and, therefore, possibly have 
higher powers of the infrared divergency. Therefore, we use scaling arguments to estimate the critical exponent for 
the ground-state energy. Below we will take care only of powers of divergencies and omit numerical factors. 

Let us assume that the matrix elements of the perturbation operator Va between low- lying states jV'n) involved into 
the PT (having = but different total 5*) at ^ oo behave as 

{4,,\Va\^j) a (20) 

Collecting the most divergent parts in all orders of the PT, the correction to the ground state energy takes a form: 

oo 
m=0 

where Cm are unknown constants and 



Ek — En 



aN^ (22) 



is a scaling parameter, which absorbs the infrared divergencies. 

The scaling function fa{x) at a; —> is given by the first order correction. In the thermodynamic limit (x — > oo) 
the behavior of fa{x) is generally unknown, but the natural condition Eg^a) ^ N at N —^ oo requires 

U{x) ^ ^ (23) 

and, finally 

E^ia) ^ -Not'l'^ (24) 

The obtained expression is in agreement with the exact equation (I13|) for the ground state energy, which justifies 
our assumption about the behavior of the matrix elements (|20p . 

Moreover, exploiting the fact that the system in the region A < 1 is in a spin-liquid phase, the correction to the 
ground-state energy has a form [3l| 

Eq = iVeo (25) 

where eo is the ground state energy at iV ^ oo and the central charge is c = 1 in our case. 

In order to reproduce such 1/A^ correction to the energy, the asymptotic of /a (a;) at large x should have a form 

/„(x)=aV^+4= (26) 
\Jx 

with some constants a and h. 
So, from Eq. ([2T|) we find 

£;o(a) = ^Nao?l^ - (27) 



and, therefore. 



^^sound ~ Va (28) 



which agrees with an exact result Wsound = \/ otj'^ at a ^ 1 32]. Thus, the scaling estimates give us the correct 
exponent for the sound velocity as well. 

Now let us consider the PT (|16p containing both channels Ya and Vj. In order to estimate the powers of divergency 
of high-orders in this PT one needs to know the A^-dependence of the matrix elements (V'i| Vj . In general, it is 
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unknown. However, one can restore these matrix elements from the known exact expression for NNN spin correlator 
in the ground state \ipo{a)) at J = and some small value of a |33j: 



or, in other words, 



(^o(a)l (S„ • S„+2 - i) \M»)) = -^^^'^ (29) 



(i^o(a)l Vj |^o(«)> = -^a3/2 JiV (30) 
37r 



On the other hand, collecting all contributions of the PT to the linear term in J similar to done in Eq. (j2ip . we 
arrive at a scaling form in small parameter a 

(Vo(a)l Vj IV'o(a)) ^ m Vj • fj{x) (31) 

with X = aN"^. 

The comparison of Eqs. ([5(I)) and PT|) immediately leads to the results 

{^^\VJ\^PJ) ^ JN-^ (32) 

and 

fAx) ^ (33) 
So, the matrix elements ("0^1 Vj \ipj) are small enough to eliminate dangerous denominators: 



Ek — Eq 



J (34) 



which, in turn, implies the absence of infrared divergencies in Vj channel. Thus, the perturbation Vj does not form 
a scaling parameter and the ground state energy has regular expansion in J. 

It is natural to expect that the behavior of the matrix elements of the type ((32|l remains the same up to the point 
J = 1/4. It results in the expression for the ground state energy at J < 1/4: 

Eo = -Na^^^gjiJ) (35) 

where gj{J) is some unknown smooth function, which at small J has the expansion in accordance with exact results 
(HI, (EOl): 

iVa3/2 

However, at approaching to the point J = 1/4 one should take into account that the excitation spectrum is 
Efe = (i — 2J) k"^ and the excitation energies become 



Ek~E,^^ (37) 



This modifies the scaling parameter 



and the expression for the energy 



(38) 



Eo (39) 



Similarly, the sound velocity at J 1/4 behaves as 



t'sound ~ \/ ^ - J^/a (40) 
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V. PERTURBATION THEORY NEAR THE TRANSITION POINT J = 1/4, A = 1 

At J — 1/4 and A = 1 the ferromagnetic ground state becomes degenerated with a singlet spiral state [l^. At 
A < 1 the ground state obviously lies in the S'^ = sector. Therefore, in order to determine the transition line 
between the phases I and III one should develop the PT both to the ferromagnetic state with S'^ = and to the 
singlet spiral state. 

A. The PT to the ferromagnetic state with = 

Let us represent the Hamiltonian in a form 

H = Ho + Vs + V^ 

Ho = - "^{^n ■ S„+l ^ |) + I X!^^" ■ ^"+2 ~ 4) 

- 7^(S„.S„+2- J) (41) 

We assume that the behavior of the matrix elements of operators Va and Vy remains the same as in the region 
J < 1/4 (see Eg. d^ni) and Eq. (|5^ ). However, the scaling parameters are modified due to the changing in one-particle 
excitation spectrum, which is Ek — at J = 1/4. So, the low- lying excited states involved in the PT ((4T|) behave 
as 

Ek -Eo- N-^ (42) 
Now according to Eqs. ip^ . ([M)) both channels and Vy produce the scaling parameters: 

X = aN^ 

y = jN^ (43) 
Thus, as follows from Eq. (PT|) exactly at J 1/4 (y = 0) the ground state energy can be written in a scaling form 

Eo{a) = -Na'^^fix) (44) 

This scaling and the critical exponent is confirmed by numerical calculations, where the function f{x) is calculated 
on finite chains with different N and a for the ground state with fc = (see Fig. 2). As one can see on Fig. 2, all data 
lie perfectly on one curve f{x) and in the thermodynamic limit the function f{x) — > 0.08. We show in Fig. 2 that the 
same scaling pi)) is valid for the lowest excited state with /c = tt as well and that the corresponding scaling function 
has the same thermodynamic limit f{x) 0.08. 

The system at J = 1/4 and a > is in a spin-liquid phase, which is verified by 1/N behavior of low-lying excitations 
calculated on finite chains (see Fig. 3). From the scaling equation we can extract also the critical exponent for 
the sound velocity: 

«sou„d ~ a'/^ (45) 
In case when both Vq. and Vy play simultaneously, the scaling estimates (|43p give 

Eoia,j) = -Na'/^f{x,y) (46) 

In the thermodynamic limit, when both a; — > 00 and y — > 00, the scaling function f(x,y) becomes a function of one 
variable (independent on N) 

2 2 

^ ^ y_ ^ (47) 

X a 

and the ground state energy takes a form 

Eoia,j)^-Na^/S{iy) (48) 
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FIG. 2: The scaling function f{x) in Ea. (|44p for the ground state energy and the lowest excited state at J = 1/4. 



B. The PT to the singlet spiral state 

The PT to the singlet spiral state with pitch angle (p in the isotropic case a = was developed in Ref. [131, where 
it was found that the energy has a scaling form 

2 

E{0,^,^)^-N^ + Nip^fi^N^ipN) (49) 

where the first term comes from the first order of the PT in 7 and the second one originates from the scaling estimates 
of the infrared divergencies of higher-orders of the PT. 

Comparison of Eq. (|^5)) and Ea. (|^5)) leads to a general expression for the energy, which correctly reduces to both 
cases at ^ and a ^ 

2 

E{a, 7, if) = ~N^- + Nip^f{aN^, -fN^,<pN) (50) 

(in fact, this equation can be derived in a similar manner as was done in Ref . [33| ) 

In the thermodynamic limit, when all variables in the scaling function in Eq. (|50p tends to infinity, the scaling 
dependence transforms to a function of two variables 



„2 



Eia,j,ip) = -N^+Nip'g{^^,iy) (51) 



where 



a 

A* = — 

v,2 



1 (52) 
a 
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FIG. 3: A^-dependence of the energy gap between the ground state and the lowest excitation of the model ([2]) at J = 1/4 and 
A = 0.96. 



Generally, the function g(/i, v) is unknown. However, we can identify some of its properties. At first, in the limit 
(y9 — > we should reproduce Eg. (j48| . Moreover, in the spin-liquid phase the spiral states with ^ N^^ should 
describe sound-like excitations with the sound velocity (|45|) . These requirements suggest that in the limit /i ^ oo the 
function g[^,v) has an asymptotic behavior 

hm .9(m, v) ^ -fi'^^giiiy) + ^^I^g2(v) + o(^^^l^) (53) 

11 — 'OO 

One can check, that this expression reproduces the sound-like excitations at 7 = 0: 

Q,3/4 

£;(a,0,^)-i?(a,0,0)^iVa3/V^ (54) 



In the limit /i ^ 0, according to Ref.[3J| we have 



lim (7(/i, v) ^ 0(1) 



(55) 



where constant A describes the excitation spectrum at the transition point a = and 7 0. Finite-size calculations 
give for this constant the value A w 0.0065. 

Summarizing all above we extract explicitly the corresponding terms and obtain the following expression: 

J-i?(a,7, (^) = - ^ - a^'I^M + o?l^p^g2(y) + V + ^'53 (m, ^) (56) 



where the function g3{fi,i') has limits 



limgsin,!^) - 0(1) 

fj,—>0 

lim gsil^,!^) ~ o(/z^/'') 

fl.—*00 



(57) 
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The minimization of the energy (|56p over the pitch angle </? gives equation for ip-^ 



dip 



= 



or after some algebra 



(58) 



(59) 



We see that in the isotropic limit a ^ (^"^ (/i ^ 0) the pitch angle is defined by the constant term in right-hand 
side of Eq.§^ (using Eq.fSTl)) 



\5AJ 



(60) 



which reproduces the result of Ref. [34| . 

In order to find the commensurate-incommensurate transition line, where the pitch angle ipmin vanish, it is more 
convenient to rewrite Eq. ([59|) in a form: 



7 



v3/4 



232 (i^) 



1 



,3/4 



5A + 553(Ai,i^) -4a* 



dfi 



(61) 



From Ea. ((6T|l and Ea. ([57| one can see that the right-hand side of Eq. (|6T|) tends to zero at /i ^ cx), which corresponds 
to the limit ip 0. The left-hand side of (|6ip is independent on ip and vanishes on the transition line: 



7 = 252(0)0^/* 



(62) 



(we note, that on the transition line ((62| — 0). 

Hence, at approaching to the transition line ([62|l the pitch angle Pmin smoothly goes to zero. So, the line ((62|) 
determines the second-order transition line between the commensurate spin-liquid phase I with p — and the 
incommensurate spiral phase III with p> ^ 0. 

Another question that can be studied concerns the low-lying excitations in the incommensurate phase. According 
to Eq. (l58|) the behavior of the energy near pmin is expanded as: 



E{a, 7, (p) = E{a, 7, </3,nin) + 



((yf - Vminf d'^E{a,'J, p) 



2 9(^2 

The second-order derivative of the energy at Lp — pmin can be estimated as 

(92-B(q,7, Pniin) 



dp"^ 



Thus, the states with 



describe gapless excitations with the energy 



27r, 



(63) 



(64) 



(65) 



(66) 



Certainly, there is no helical LRO in the spiral phase and the spin correlations decay on large distances. However, 
the nature of the spiral phase manifest itself in the incommensurate position ^max of the maximum of structure factor 



S{q) = • 



(67) 



When the 0(3) rotation symmetry is broken by the anisotropic term Va-, the incommensurate nature of the spiral 
phase remains in the x-y plane. So, in this case we associate the pitch angle of the spiral pi with the position of 
maximum of the structure factor ^max in the easy plane 



(68) 
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The numerical calculations on finite chains show that for a fixed small value of 7, gmax decreases via consecutive 
sharp jumps on the value ^ from some finite value at a = to zero on the transition line 

a « 13.97^/3 (69) 

where the incommensurate phase III terminates (see Fig.l) and the transition into commensurate spin- liquid phase 
takes place. Thus, the numerical calculation confirms the found critical exponent for the transition line (|62p . The fact 
that the pitch angle ip tends to zero at approaching to the transition line ensures that this line is the second-order 
transition. 



VI. THE EASY-AXIS CASE 



In the easy axis case for J < 1/4, the fully polarized state ||| . . . |) is evidently the ground state. In the region 
J > 1/4 one should compare the fully polarized state energy with the energy of the spiral state. The finite-size 
numerical calculations show that for a fixed small 7 the increasing of easy-axis anisotropy leads to the decrease of the 
pitch angle but the ground state remains in the sector with total — 0. At some critical value of Ac the transition 
from the state with S*^ = and some finite value of the pitch angle to the fully polarized state occurs. Thus, in 
contrast to the easy-plane case, the transition from the spiral phase to the ferromagnetic phase is the first-order one. 

The finite-size numerical calculations also show that in the spiral region 1 < A < Ac, it is sufficient to take into 
account only the first-order correction in (A — 1) to the spiral state. That is the energy of the spiral state is 

E,p - -aN^''/^ - (70) 

and the transition to the fully polarized state with the energy 

Ef = (71) 

takes place at 

Ac = 1 + 607^/3 (72) 

Unfortunately, the finite-size calculations do not allowed to find the factor a in Eq. ([7^ because of irregular behavior 
of Ac with N . However, we believe that the MFA gives a good estimate for this transition line pT|). 



VII. SUMMARY 



We have studied spin-1/2 zigzag chain with weakly anisotropic ferromagnetic nearest-neighbor and antiferromagnetic 
next-nearest-neighbor interactions. It was shown that the ground state phase diagram consists of three phases: the 
fully polarized ferromagnetic phase, the commensurate spin-liquid phase and the incommensurate phase. Thus, the 
incommensurate phase established for the isotropic case survives weak anisotropy of interactions, though in this case 
the incommensurate nature of the ground state reveals itself in the x — y plane. 

Using scaling estimates of the infrared divergencies in the perturbation theory we obtained the scaling expression for 
the ground state energy both for commensurate and incommensurate phases. This allowed us to determine non-trivial 
critical exponents in the behavior of the phase transition lines, which were confirmed by finite-size calculations. We 
found also that in the easy-plane case the transition from the commensurate spin-liquid to the incommensurate phase 
is of the second order one, while in the easy-axis case the transition from the fully polarized state with 5^ — 5'f„^^ to 
the incommensurate state with 5^ = is evidently of the first order. 

In this paper we have focused on studying of the model ([2]), which is a particular case of more general model ([T]). 
However, the obtained results for the model (O remain valid at least qualitatively for the model H]). As an example 
we present in Fig. 4 the phase diagram near the transition point J = 1/4 of the model ^ in the case Ai = A2. We 
see that the phase diagram in this case is very similar even quantitatively to that shown in Fig.l. 
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